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Extraction of timelike entanglement from the quantum vacuum 
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Recently, it has been shown that the massless quantum vacuum state contains entanglement be- 
tween timelike separated regions of spacetime, in addition to the entanglement between the spacelike 
separated regions usually considered. Here, we show that timelike entanglement can be extracted 
from the Minkowski vacuum and converted into ordinary entanglement between two inertial, two- 
state detectors at the same spatial location — one coupled to the field in the past and the other 
coupled to the field in the future. The procedure used here demonstrates a clear time correlation 
as a requirement for extraction, e.g. if the past detector was active at a quarter to 12:00, then the 
future detector must wait to become active at precisely a quarter past 12:00 in order to achieve 
entanglement. 

PACS numbers: 03.70.+k, 03.65.Ud 



I. INTRODUCTION 



The quantum vacuum is theorized to exhibit a variety of thermal effects associated with spacetime horizons, in- 
cluding Hawking radiation in black hole spacetimes [T], Gibbons-Hawking radiation in de Sitter spacetime [5], and 
Unruh-Davies radiation for accelerated observers in Minkowski spacetime [31 3] . The predicted thermalization of the 
vacuum in these examples is associated with intrinsic non-classical correlations, i.e. quantum entanglement, between 
different regions of spacetime. The accessibility of Minkowski spacetime has given rise to research interest in the 
potential for direct manipulation of this field entanglement [5HS]. However, given that the indirect signature of the 
entanglement, i.e. the thermalization of the vacuum, has remained well outside the observable regime, the prospects 
^ for direct observation of the underlying entanglement have seemed remote. 

Very recently, we have predicted a new thermal effect closely related to that of Unruh-Davies, which affects an 
inertial particle detector switched on at i 0, having a detector energy gap which is scaled with time as ^ , where a 
^— I is the scaling constant 10 . Such a detector registers a thermal response to the vacuum, with temperature T ~ 
^ The requirements for observing this effect appear to be much closer to current technological ability than for the 
Unruh-Davies effect which requires very high accelerations. The origin of the thermalization in this case is that the 
Minkowski vacuum of massless quantum fields is in fact entangled between timelike separated regions of spacetime, 
in close analogy to the case of spacelike separation usually considered [TT] . 

Here we demonstrate the practicality of extracting this entanglement and hence of directly observing the entangle- 
JT^ ment of the quantum vacuum for the first time. The conceptual novelty of quantum non-separability across time raises 
^— ^ some immediate questions about its relation to standard entanglement, which is usually imagined to be a property 
^ of a quantum state at a particular time. We explore the relation to ordinary entanglement by demonstrating that 
• . timelike vacuum entanglement can be extracted and converted into ordinary, constant-time entanglement between 
. ^ detectors, one which couples to the field in the past region P, and another which couples to the field in the future 

[v>( region F of Minkowski spacetime (see Figure 1). We find that the timelike nature of the entanglement leads to peculiar 
; I correlations in time. 

In section 2, we review the specific state of the field in F and P, noting entanglement and the thermal state of the 
field when restricted to F or P alone. We also review the thermal single-detector response of the ^ scaled detector. 
In section 3, we describe the two-detector entanglement extraction procedure, and a description of the conditions 
under which entanglement extraction is possible — the central result of this paper. Section 4 contains our concluding 
remarks. 
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FIG. 1: Spacetime divided into quadrants consisting of regions contained by the future and past light cones (F and P), and 
the right and left Rindler wedges (R and L). 



II. TIMELIKE NON-SEPARABILITY IN THE MINKOWSKI VACUUM STATE 



Before discussing the procedure for extracting timelike entanglement from the vacuum, we first recall the form 
of the Minkowski vacuum restricted to F-P, and see the entangled state explicitly, for the simplified case of a 1+1 
dimensional spacetime. We examine the left moving sector of a free massless scalar field, for which the left and right 
moving sectors may be quantized independently, and we use the approximation of a discrete set of modes labeled by 

LUi. 

For entanglement to be defined, we require the field in the future to be quantized as an independent system from 
the field in the past. For a scalar field (j){x), this is satisfied when [(f){x f) , (t>{x p)] = 0, where xp and xp are any 
timelike separated points. This condition is satisfied for massless fields, and in the limit of an arbitrarily small but 
non-zero mass the commutator is bounded by a maximum value of ^ for any timelike separated points [121 113| , and 
thus a sufficiently small mass will also allow an approximation of independent systems between F and P. 

It is well-known that in 2-d spacetime, the Minkowski vaccum restricted to the right, R, and left, L, Rindler wedges 
(see Fig. 1) can be expressed as an entangled state of the Rindler modes in the following way, which can be understood 
as the basis for the Unruh effect [T^ : 

|0m> = E ^, («5f«f;!)"'|0fl). (1) 

i ni=0 

where \Qr) is the Rindler vacuum, and a^'^ is the creation operator for a right Rindler particle, corresponding to the 
solution g^{T + e) = (47rti;)^^/^e~"^'^'^'^^^ in the coordinate system t ~ a^^e"'^ sinh(ar), z = a^^e"*^ cosh(ar), with t 
and z the usual Minkowski coordinates (and analogously for af^^). 

The Minkowski vacuum restricted to F-P takes an exactly symmetrical form, when expressed in terms of the 
"conformal modes" and [TU] : 



iom) =n^^ E ■ 

i iii—O 



(«5;«r/)-|0T). 
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(2) 



Here, F is coordinatized \>yt — a ^e"^^ cosh(aC) and z — a ^e'"' sinh(aC) with P coordinatized hyt = —a ^e"'' cosh(aC) 
and z = —a^^e°-^smh{aQ. The creation operators correspond to the solutions g^{r] + C) = (47ra;)^^/^e^'"(''+''^ and 
gP{fj + C)^ (47rw)-i/2e-*'^(^+f\ with |0t) their vacuum. 

This symmetry reflects the fact that the conformal modes in g^ are in fact the same solutions as g^, continued 
from R into F. 

Also mirroring the case of spacelike entanglement, the state of the field in F (or P) alone is a thermal state of the 
conformal modes: 



jii=0 



2TrniUJi/a\ F 



(3) 
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In 3+1 dimensions, an inertial Unruh-DeWitt detector (that is, one moving along the trajectory x = y = z = 0, 
t = a~^e'"') can also be seen to respond to the vacuum in a manner identical to that of an accelerating detector, 
provided it is designed to evolve in the conformal time 77, rather than in the proper time r. This means that the free 
Schrodinger equation reads: 

= Ho^ (4) 

where the eigenvalues of Hq are taken to have a constant gap, E. Transforming this equation to Minkowski time t 
reads: 

4v|/=^vi,. (5) 

dt at ^ ' 

In other words, we are describing a detector whose energy gap must be scaled as ^ in ordinary Minkowski time t. 

The interaction term Hj may be taken to be the standard Unruh-DeWitt term, which for a two-state detector takes 
the form Hj = a4>ix{t))[\0){l\ + |1)(0|]. We thus consider the fuU Hamiltonian: 

i|^* = (Ho + e"''^^/)* (6) 

The interaction term acquires the exponential factor due to the change of variables to conformal time, and because 
the coupling to the field is the standard one, and not scaled in time. 
From here, one finds the detector reponse function to be: 

/•CO poo 

F{E)= dT] d7/'e-'-^(''-''')e''(''+''')i:>+(77,r;') (7) 



where D+{7j,r,') = {OM\Hv)Hr]')\OM)- 

The limits of integration correspond to a detector which becomes active at t — 0. Another symmetry between the F- 
P case and the R-L case now becomes important — the two-point function along the inertial trajectory x — y = z = 0, 
t = a^^e"'' can be calculated to take the form: 

2 -a(jj+r)') 

D+ (77, ry') = ^ (8) 

^' " 4sinh2(|(77-,7')-«e) 

while the two-point function along the trajectory of a uniformly accelerated trajectory t = sinh(aT), x — y = 0, 
z = cosh(ar) takes the form: 

D+ (t. t') = . . (9) 

^ ' ' 4sinh2(|(r-T') -ie) ^ ' 

This symmetry leads to a formally identical response function integral in the two cases (the additional exponential 
factor canceling with an exponential factor in the response function integral in the F-P case), and thus to the same 
thermal signature when the integral is evaluated by standard techniques IT5] . In the case of an accelerated 
trajectory, the acceleration plays the role of temperature, giving T(j = 5^^, with one degree Kelvin corresponding to 
an acceleration on the order of 10^" p^. In the inertial case, it is the magnitude of the scaling of the detector energy 
gap, a, which plays the role of temperature, giving T — with one degree Kelvin corresponding to a scaling on 
the order of 100 Gigahertz. In the following section, we will consider two such scaled detectors — one in the future 
and one in the past. 



III. TIMELIKE ENTANGLEMENT EXTRACTION 



Here we present our main findings, in three components. First, we describe the extraction procedure in some 
detail, and show that extraction of timelike entanglement is possible. Second, we describe a basic symmetry property 
of our procedure, namely that identical entanglement can in principle be extracted between regions of arbitrarily 
great timelike separation, but that larger timelike separations require a correspondingly longer interaction time for 
the detectors to achieve the same degree of entanglement. Third, we show a fundamental time correlation in the 
extraction procedure. 
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A. Entanglement Extraction Procedure 



We now consider two two-state, energy scaled Unruh-DcWitt dcttectors, one of which is active in F, with the other 
active in P. Due to the properties of the field commutator (Pauli- Jordan function), for which massless fields satisfies 
[(j)(x) , 4>iy)] = for timelike separated points x and y, the detectors thus interact entirely with independent systems. 

In F, the detector moves along the inertial trajectory parameterized by conformal time rieLSx = y = z = 0, 
t = a~^e'^'^, while in the past, P, the trajectory is parameterized hyx = y — z — 0, t — —a~^e'^^. The detectors will 
be sensitive to the frequency E with respect to the conformal time in their respective quadrants. This requires the 
following energy scaling in terms of the Minkowski time t: 

Hf - ^+Hi (10) 
at 

Hp = ^+Hj (11) 
at 

where Ho\0) — Eo\0) and Ho\l) = Ei\l), with Ei — Eq = E, which is taken to be positive. The minus sign appearing 
in Hp cancels with the negative value of t in the past, so that the interpretation of "ground state" |0) and "excited 
state" |1) is standard for both detectors. 

At t — — oo, we take the state of the detectors to be |00), and we wish to determine the state at i = oo — 
specifically, we will be interested to know whether the state of the detectors is entangled. Given that the energy 
gap of the P-detector diverges as it approaches i = 0, we make the following assumption: after the P-detector has 
interacted with the field in P, but before t — 0, we assume that the detector energy scaling is adiabatically turned 
off (rather than allowed to "blow up" at i = 0), so that the state of the P-detector is effectively "frozen" after its 
interaction with the field. Similarly, the energy scaling of the future detector is adiabatically turned on prior to its 
interaction with the field. 

To determine entanglement, we follow an approach analogous to that of Reznik, Retzker, and Silman IL, who 
first studied vacuum entanglement extraction from the spacelike separated separated Rindler wedges using two, two- 
state detectors. Specifically, we look for a positive value of the negativity of the two-detector state at t = oo, which 
is the necessary and sufficient condition for the non-separability of the 2x2 dimensional system formed by our 
detectors ,16. ,17) . 

To express the state at t = oo, we use perturbation theory in the conformal time, the same as in the single-detector 
case, but we now include the "window functions" XF{fl) s-i^d xp{f])i defined in F and P respectively, describing the 
interval over which the detectors are active. 

To second order, the state at i = oo thus takes the following form: 

I*) = (1-C)|0m)|00) 

-z r d77XF(»?)e'^''e-*^V(^)|OA/)|01) 



d77Xp('7)e'^^e^^V('7)|OA./)|10) 

/oo 
d^XF(??)xp(??)e"(''+'^^e-'^("-'')0(ry)0(r7)|OM)|ll) (12) 
OO 

To simplify notation, we define the following (unnormalized) field states: 

/oo 
dr;xF(?7)e'^''e-'^V(^)|0M) (13) 
-OO 

/oo 
df]xp{fl)e''^e^'''$mOM) (14) 
-OO 

/oo nOO 
dv / rff?XF(r?)xp(??)e'^(''+'')e-'^("-'5)0(,7)0(r7)|OM) (15) 
-OO J — OO 

The state at t = cx) can then be written in the following simplified form: 



1*) = (1 - C)|Om)|00) - i\AF)\01) - i\Ap)\m) - |X)|11) 



(16) 
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We now trace over the field degrees of freedom, and obtain the two-detector density matrix in the basis |00), |01), 
|10), 111). 



( N -{X\{)m) 

{Af\Af) -{Ap\Af) 

-{Af\Ap) {Ap\Ap) 

V-(OmI^) {X\X) 



(17) 



where N = 1- {X\X) - {Ap\Ap) - {Af\Af). 

To lowest nontrivial order, the negativity M{p) is given by: 



N{p) = |(Oa/|X)| - ^{Af\Af){Ap\Ap) 



(18) 



The condition for non-separability of the two-detector state is that N{p) > 0. When the window functions xf 
and xp ^^'^ symmetrical about t = (as well as the scaling constant a and conformal frequency gap E, so that 
{Af\Af) = {Ap\Ap)), the non-separability condition amounts to the following: 



> 



dv / d??XF(r/)xp(^)e'^(''+'')e-^^("-'?)(OM|<^(??)<^(^)|OM) 
dv / rf^'Xi=^(^)Xi=^(^')e°^''+'''^e-'^("-''')(OM|<^(?7)0(?7')|OM) 



(19) 



We have used the notation that (p{r]) = (f){x{ri)) along the trajectory t — a ^e"''',x = in F, while in P wc have 
that t = -a^^e^^f = 0. 

The quantity {0m\4>{x)(I){x')\0m) takes the ordinary regularized form — [{t — t' — ie)^ ~ ~ x')"^] [T^, and 
thus for the given inertial trajectories, a coordinate transformation yields (for appropriately rescaled infinitesimal 
regulator e): 



{OM\mkv')\OM) = 



167r2 sinh2( 



2( a{r)-ri') 
_^2g-a(r,+f;) 



IGvr^ cosh (- 2 
We now consider a specific paired set of window functions: 

Xf{vi) = e"''' 

-2 

The entanglement condition A/" > thus reduces to the following: 



drj 



> 



drj 



dr7e-"'-*''e-^^(''-*'' cosh" 



(20) 
(21) 



(22) 
(23) 



(24) 



We express this compactly as Ix > Ia- The first quantity, Jx, is not singular, and may be numerically integrated 
in straightforward fashion with Mathematica. For E = 1 and a = 2, the value of Ix is approximately 1.561. The 
second quantity, I a, is formally identical to the (unnormalized) response function integral of an accelerated detector 
with a window function, where E would represent a fixed proper-energy gap, and the integration variable 77 would 
represent the proper time (rather than the conformal time it represents here). 

To make I a convenient to compute, we first make use of the identity csc^(7ra:) — X]feL-oo {x-k)^ ■> ^'^ that we can 
write: 
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/oo poo 



' — OO <J — OO 

poo poo 



OO J — OO 



(a(?7 — 77') — ie)^ ^—J (a(77 — ?/') + ink)^ {a{r] — r/') — inky 



00 /'OO 



/incrtial + / d,?' e""^ - e-^^^(''-'/) ^ f + ^ f . (25) 

■/-co \^ (a('7 - ) + «7rfc)^ (a(?7 - r;') - ^tt/c)^ ^ 

The sum/integral in the final line can be evaluated numerically in Mathematica, while the term /inortiai has been 
evaluated by Satz |18j with careful attention paid to regularization, since it is formally identical to the response 
function of an inertial detector with fixed energy gap E. For arbitrary window function, it takes the form: 



/OO /"OO /"OO 

-00 J —00 Jo 

roo i„ ^00 



1 COs(i?sJ 



S 



2/ ^/ dr]xiv)[xiv)-xiv-s)] (26) 



* J-oo 



where a change of variables was made to s = 77' — rj" . Numerically evaluating this quantity in Mathematica, and 
combining with the sum/integral term for our chosen window function xpiil) = sT^ yields a value of I a which is 
approximately 1.273 (i.e. smaller than Ix ~ 1.561), and thus yields a positive value for the negativity, demonstrating 
the nonseparability of the detectors at t = 00. 



B. Time Translation of the Window Functions 



The value of the negativity calculated in the previous manner is invariant under a simultaneous translation of the 
window functions in the conformal times rj and 77 by an amount x to the window functions xpiv) — e^^"^^^^ and 
Xpiv) = e"'-''"^-' . However, different values for x result in different widths of the window functions in the Minkowski 
time t (see Figure 2). As the window functions are shifted away from the t = origin, the detectors must remain 
active for far longer, in order to achieve the same degree of entanglement. This supports the intuition that most of 
the field entanglement is concentrated in the region close to the edges of the light cone. 



X 




FIG. 2: Three symmetrically paired sets of window functions, Xf{v) = e~^^~^'' and Xpiv) = e"'*'"^-' , plotted in Minkowski 
time t, which produce identical entanglement between the detectors. As the window functions are shifted away from t = 0, the 
detector must remain active for far longer to achieve the same degree of entanglement. 



7 



1.0 r 



0.5 - 




-1.5 - 



-2.0^ 

FIG. 3: Numerical value of the quantity Ix — Ia, for choice of window functions xp = e ^ and xf — e ^' , for detector 
parameters E = 1 and a = 2. The detectors are entangled at f = oo for positive values of Ix — I a only. The position of the 
window functions must by sufficiently symmetrical about t = (corresponding to a sufficiently small value of x) in order to 
extract entanglement from the vacuum between the timelike separated regions F and P. 



Quantitatively, the price to be paid in terms of the total "volume" of interaction time to entangle the detectors 
to the same degree (as measured by Minkowski time) is exponential in the dispacement x in conformal time. That 
is, the ratio XFiv ~ x)dt/ XF{''i)dt — e"^. This is analogous to the exponential fall-ofF of entanglement with 
distance in the spacelike separated case which was observed by Reznik, Retzker, and Silman [11) . 



C. Time Correlation 



The negativity does not remain constant, however, if only one of the window functions is shifted, spoiling the 
symmetry about t = Q. In fact, the entanglement can be completely destroyed by shifting one of the window functions 
sufficiently far away from its symmetrical location in time. One can readily verify that the quantities {Ap\Ap) and 
{Ap\Ap) in equation 18 are each independently invariant under a translation XF(ji) — > xf(ji — x) and X-f(^) ^ 
Xf{j] — x), and thus the quantity I a in equation 24 remains unchanged under the translation of a single window 
function. However, the quantity |(0m|^)| in equation 18 does not possess this symmetry if only one of the two window 

functions are shifted, and thus the resulting value of Ix = 



dv Coo dv XFiv - a:)xp(^-)e-^(''-'') cosh"^ (^) 



will determine whether or not entanglement has been extracted. We illustrate this by plotting the quantity Ix — Ia, 
which is proportional to the negativity. The detector state at t = oo is thus separable only when Ix — Ia is positive. 
The window functions plotted are the following: 

XF = e-^"-")' 

XP = e-^\ (27) 

This corresponds to a situation where the past detector has already interacted with the field at time ^ = 0, and 
we must now select the time at which the future detector will be active, with the quantity x signifying how far (in 
conformal time) we move away from the point of symmetry. Ix — Ia is plotted as a function of the choice of x in 
Figure 3 — positive values correspond to an entangled final state, while negative values correspond to a separable 
final state. Clearly, entanglement is maximized around the symmetrical point in time (corresponding to a; = 0), while 
a sufficiently non-symmetrical choice for x can kill the extraction of timelike entanglement entirely. 

Stated more dramatically, a detector which is switched on and off in the vicinity of a quarter to 12:00 can become 
entangled with a detector interacting with the field at the same spatial location in the future, but only if the later 
detector waits to be switched on and off at a quarter past 12:00. 
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IV. CONCLUSIONS 



We have attempted to partially answer the question of how timelike entanglement in the Minkowski vacuum can 
be related to the familiar "entanglement at a given time" between simple, two-state systems by showing that timelike 
entanglement may be extracted and converted into ordinary entanglement between two two-state detectors. We thus 
conclude that timelike entanglement may be regarded as a non-classical resource in a manner analogous to the spacelike 
entanglement that is often studied in the Minkowski vacuum, since any quantum information theoretic protocol may 
utilize conversion of timelike entanglement to spacelike entanglement as a step in the protocol. 

As a thought experiment to illustrate this possibility, we imagine a quantum teleportation protocol in which the 
entanglement resource is between a detector interacting in P, and a detector interacting in F, and all operations 
on the P-detector and the qubit-to-be-teleported take place before t — 0. Classical information alone is then sent 
into F, where the F-detector must interact with the field there at a particular time to form the other half of the 
entanglement resource. The classical information from P is then used to transform the F-detector into the teleported 
qubit. Such a protocol might be called "teleportation in time," since there exists a period after t = but before the 
future interaction time where it is not possible to recover the teleported qubit. 

We thank Nicolas Menicucci for stimulating discussions. In addition, we thank the Defence Science and Technology 
Organisation (DSTO) for their support. 
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